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ABSTRACT 

We present observational evidence of environmental effects on the formation and evolution of planetary systems. We combine cata- 
logues of resolved protoplanetary discs (PPDs) and young stellar objects in the solar neighbourhood to analyse the PPD size distri- 
bution as a function of ambient stellar density. By running Kolmogorov-Smirnov tests between the PPD radii at different densities, 
we find empirical evidence, at the > 97% confidence level, for a change in the PPD radius distribution at ambient stellar densities 
£ 10 3 ' 5 pc -2 . This coincides with a simple theoretical estimate for the truncation of PPDs or planetary systems by dynamical 
encounters. If this agreement is causal, the ongoing disruption of PPDs and planetary systems limits the possible existence of planets 
in the habitable zone, with shorter lifetimes at higher host stellar masses and ambient densities. Therefore, habitable planets are not 
likely to be present in long-lived stellar clusters, and may have been ejected altogether to form a population of unbound, free-floating 
planets. We conclude that, while highly suggestive, our results should be verified through other methods. Our simple model shows 
that truncations should lead to a measurable depletion of the PPD mass function that can be detected with ALMA observations of the 
densest nearby and young clusters. 

Key words. Planets and satellites: formation - Protoplanetary disks - (Stars:) circumstellar matter - Stars: kinematics and dynamics 
- (Galaxy:) open clusters and associations: general 



1. Introduction 

For the past decade, exoplanetary systems are being discovered 
at a spectacular rate (e.g. Mayor et al. 2004; Borucki et al. 201 1). 
Driven by these discoveries, there is an increasing interest in the 
global properties of planetary systems, from the epoch of their 
formation in protoplanetary discs (PPDs) to their long-term sta- 
bility. While there is a natural focus on internal processes that 
govern the evolution of such systems (e.g. Lee & Peale 2003; 
Dullemond & Dominik 2005; Gorti et al. 2009; Morbidelli et al. 
2009; Brasser et al. 2009; Blum 2010; Williams & Cieza 201 1), 
it is clear that not all planetary systems form in isolation and 
environmental effects should be considered as well. Particularly, 
theoretical studies show that external photoevaporation (Scally 
& Clarke 2001; Adams et al. 2004, 2006; Fatuzzo & Adams 
2008) and dynamical interactions (Bonnell et al. 2001; Pfalzner 
et al. 2005; Olczak et al. 2006, 2010; Spurzem et al. 2009; 
Lestrade et al. 201 1; Dukes & Krumholz 2012; Parker & Quanz 
2012; Bate 2012) can lead to the truncation of PPDs and plane- 
tary systems. 

While the external photoevaporation of PPDs has been stud- 
ied observationally (O'dell et al. 1993; Robberto et al. 2008; 
Rigliaco et al. 2009; Mann & Williams 2010), there is no con- 
clusive evidence of dynamical effects (Eisner & Carpenter 2006; 
Olczak et al. 2008; Reche et al. 2009). In part, this is likely 
due to the relatively short lifetimes of PPDs (up to ~ 8 Myr, 
e.g. Haisch et al. 2001; Hernandez et al. 2008; Ercolano et al. 
201 1; Smith & Jeffries 2012) compared to the time it takes stel- 
lar encounters to have an observable effect on the disc population 



(~ 0.1-1 Gyr, e.g. Adams 2010). Previous observational studies 
on the Orion Nebula Cluster (ONC) (e.g. Eisner & Carpenter 
2006; Olczak et al. 2008) did aim to find traces of dynamical 
interactions in the population of PPDs, but lacked a sufficient 
number of sources and/or suffered from uncertainties on the disc 
mass measurements. 

We address the problem statistically by considering the sizes 
of PPDs as a function of their ambient stellar density, using sam- 
ples of PPDs and young stellar objects (YSOs) from the latest 
infrared surveys. If stellar encounters truncate PPDs by tidally 
stripping the outskirts of the discs (e.g. Clarke & Pringle 1993; 
Heller 1995; Hall et al. 1996), this should be observable above 
some characteristic ambient stellar density (~ 10 3 pc~ 3 Adams 
2010), because the encounter rate increases with density (e.g. 
Binney & Tremaine 1987, Eq. 7-61). In this Letter, we find 
model-independent evidence of a change in the PPD size dis- 
tribution for ambient stellar surface densities E ^ 10 3 5 pc~ 2 at 
the > 97% confidence level. 

2. Protoplanetary discs and their environment 

2. 1. Data selection 

To verify whether a relationship exists between the sizes of PPDs 
and their ambient stellar surface density, E, we combine exist- 
ing catalogues of PPDs and young stellar objects in star-forming 
regions (SFRs) of the solar neighbourhood. 

The data for the PPDs is taken from circumstellardisks.org 
(Karl Stapelfeldt, NASA/JPL). This catalogue gathers resolved 
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PPDs that have been confirmed and described in the literature. 
If a PPD is resolved in different wavelengths (probing different 
parts of the disc, see e.g. Lada et al. 2006), the catalogue lists 
the largest measured diameter, implying that the disc radii used 
in this work are lower limits. About 75% of the PPD radii in our 
final sample (see below) are measured at wavelengths around 
1 /xm, with only ~ 25% of the PPDs (all in low-density regions) 
being observed at mm wavelengths (see Appendix A). The cat- 
alogue contains an estimate for how well-resolved each disc is 
by listing the number of diffraction-limited beams that fit within 
its diameter. We only consider the discs for which this value is 
greater than unity. This provides us with 133 PPD sources from 
which we exclude those whose host star is classified as a main 
sequence star, weak-line T-Tauri star, or Class YSO. Sources 
at distances > 500 pc (which covers all our YSOs) are also ex- 
cluded. These criteria reduce the sample to a total of 101 sources. 

To estimate the local ambient surface density of each PPD 
source, we use publicly available near-infrared data of nearby 
SFRs (see Table A. 2 in Appendix A). The ambient surface den- 
sity of stars around each PPD is estimated as in Casertano & 
Hut (1985) - see Appendix A for a detailed explanation. The 
thus-obtained angular ambient surface densities are converted to 
physical ambient stellar surface densities, E, using the distances 
listed in Table A.2. In cases where the listed PPD distance differs 
from the distance to the region that it is a member of, we adjust 
the PPD distance and size. Since a low surface density can be 
due to incomplete YSO coverage, any discs with E < 0.1 pc~ 2 
are omitted from our analysis. Moreover, the minimum PPD ra- 
dius that can be resolved increases with distance and hence in- 
troduces a distance-dependent selection bias. To avoid this, we 
exclude any PPDs smaller than the smallest radius (~ 50 AU) 
that is resolved in the most distant region of our sample, which 
is the Orion Nebula Cluster (ONC). The final sample thus con- 
tains 67 sources (see Table A.l in Appendix A). Completeness 
does not affect the densities because the surveys of Table A. 2 are 
complete down to the hydrogen-burning limit. 

2.2. A simple theoretical estimate for the truncation of 
discs 

To interpret the data, we include a rough theoretical estimate for 
the expected truncation radii of PPDs as a function of E. This 
is obtained by combining the truncation induced by each partic- 
ular encounter with the stellar encounter rate. The derivation is 
presented in detail in Appendix B. 

We use the numerical simulations of disc perturbations in 
clustered environments by Olczak et al. (2006) to obtain the disc 
radius as a function of the encounter parameters. We convert 
their expressions for the disc mass loss to a radial truncation 
assuming that it occurs by stripping the outer disc layers to the 
Lagrangian point between both stars, and adopting a PPD sur- 
face density profile E d oc r -1 (Olczak et al. 2006). Under these 
assumptions, we write for the upper limit to the disc radius 

r d (r p ,mi,m 2 ) = J= — -, (1) 

^m 2 /m 1 + 1 

where r p is the pericentre radius at which the perturber passes, 
mi the mass of the perturbed system, and m 2 the mass of the 
perturber. This approximation is validated in Appendix B. 

The encounter radius r p is obtained from the impact parame- 
ter b, encounter velocity v, and masses m\ and m 2 by accounting 
for gravitational focussing (see Appendix B). The masses are as- 
sumed to follow a Salpeter (1955) type initial mass function in 
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Fig. 1. Top panel: Protoplanetary disc (PPD) radius distribution versus 
ambient stellar density. Squares, triangles, stars, and crosses correspond 
to ONC sources, T-Tauri stars, Herbig stars, and YSOs, respectively. 
Dotted lines represent the theoretical truncation for ages of 0.3 (top) 
and 1 Myr (bottom). Bottom panel: KS test pKS-values for the PPD 
radii at that density to be consistent with the sample at lower densi- 
ties, using {6, 7, 8} (triangles, squares, diamonds) objects per density 
bin (see text). The grey thick line shows the p K s-value for ONC PPDs 
only when dividing the subsample in two at the indicated density. The 
horizontal dashed line and grey area indicate the median p K s and its la 
dispersion retrieved from 30,000 Monte Carlo experiments to verify the 
significance of the results (see text). Dotted lines mark 95% (top) and 
99% (bottom) confidence levels. 

the range 0.1-100 M Q , and the velocity distribution is taken to 
be Maxwellian with a velocity dispersion of a = 2 km s~\ as 
is typical for SFRs (Hillenbrand & Hartmann 1998; Covey et al. 
2006). This enables the derivation of the encounter rate as a func- 
tion of b and v (Binney & Tremaine 1987), which for a given age 
provides the total number of encounters n. Because encounters 
with pericentres at inclination angles 9 il 45° with respect to 
the disc plane affect the disc only mildly, about 30% of the n en- 
counters lead to the disc truncation described by Eq. 1 (Pfalzner 
et al. 2005). We use the probability distribution functions (PDFs) 
for b, v, and m 2 to calculate the PDF of the 'most disruptive' en- 
counter, i.e. the parameter set that gives the smallest disc size, 
according to Maschberger & Clarke (2008, Eq. A5). We then in- 
tegrate the product of the disc radius r d , the PDF of the 'most 
disruptive encounter', and the mass PDF of the perturbed object 
mi to obtain the expected disc truncation radius r tr . To com- 
pare the theoretical estimate to the observations, we relate the 
stellar volume density v to the surface density E as v = E/2i?, 
where R w 2 pc is a typical radius for the SFRs in our sample 
(Hillenbrand & Hartmann 1998; Evans et al. 2009). 

3. Results 

3. 1. Evidence for environmental effects 

The upper panel of Fig. 1 shows the observed PPD radii versus 
E. The distribution is relatively insensitive to the ambient den- 
sity until E ~ 10 3 5 pc~ 2 , where it appears truncated at large 
radii. This is consistent with the simple theoretical approxima- 
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tion from Sect. 2.2, which predicts a truncation at these densi- 
ties for ages between 0.3 and 1 Myr (see Fig. 1). The affected 
PPDs are all in the ONC. If the truncation is interpreted as being 
due to dynamical interactions, the theoretical curves suggest that 
youngest sources in the ONC have ages of ~ 0.6 Myr, which is 
reasonably consistent with observations (Palla & Stahler 1999; 
Da Rio et al. 2010; Jeffries et al. 2011). However, this is not 
a unique explanation, since the truncation might also be due to 
external photoevaporation by nearby massive stars (e.g. Clarke 
2007). It is also important to keep in mind that some fraction of 
the PPDs at this projected surface density will actually reside in 
a region of lower volume density, either behind or in front of the 
high-density core of the ONC. The distribution also shows some 
evidence of a density-independent upper limit to the PPD radius 
of ~ 10 3 AU, which could be related to binarity (Artymowicz & 
Lubow 1994; Kraus et al. 2012) or be intrinsic (Basu 1998). 

To test the statistical significance of the change in the ra- 
dius distribution, we perform a Kolmogrov-Smirnov (KS) test, 
which is theory-independent and thus insensitive to model as- 
sumptions. It gives the probability p^s that two samples were 
drawn from the same parent distribution. Starting at the high- 
density end, we first divide the sample in density bins of {6, 7, 8} 
objects per bin, which represents a balance between good statis- 
tics and enough bins to resolve the regime where the radius 
distribution changes. The KS test is then carried out compar- 
ing radii in a bin at density Sbin with those at lower densi- 
ties (i.e. 10 _1 pc~ 2 < E < Sbin)- We do not include bins at 
Sbin < 10 1 pc~ 2 to avoid low-number statistics in the reference 
sample. The bottom panel of Fig. 1 shows the results of the test. 
At intermediate densities, all KS tests give high j?KS- va l ues > but 
for densities higher than E <; 10 3 5 pc -2 there is a pronounced 
drop. As such, the KS test yields a detection of a change in the 
PPD radius distribution at the > 97% confidence level. Note that 
this does not change when excluding the largest PPD in the sam- 
ple (at S ~ 200 pc" 2 ). The result also holds within the ONC 
only (grey line in Fig. 1), when dividing the ONC subsample 
in two at each density and running a KS test for the radii at both 
sides of the separation. To check the result, we performed 30,000 
Monte Carlo experiments in which the KS test was applied in the 
same way to distributions of randomly paired radii and surface 
densities, i.e. erasing any possible correlation. Figure 1 indicates 
the resulting median pxs an d its ler dispersion, showing that the 
results obtained for the original sample are unlikely to be due to 
the adopted statistical method. 

3.2. Implications for habitable zone occupancy lifetimes 

If the drop of the pKS _va l ue at high densities is indeed caused 
by dynamical truncation, then the model can be used to give the 
maximum time during which the habitable zone (HZ) can host 
a planet or a PPD ('HZ lifetime'). We calculate r tr as a func- 
tion of E and t as in Fig. 1, but without averaging over the host 
stellar mass to retain the mass dependence. We then compare it 
to the inner radius of the HZ, rjf z , and determine until which 
age they overlap, as a function of mass and ambient density. The 
radius r^ z depends on the required radiative equilibrium tem- 
perature T cq , on the stellar luminosity L and on the properties 
of the planet (i.e. /, a proxy for atmospheric thermal circulation, 
and the Bond albedo A, see Kasting et al. 1993; Tarter et al. 
2007): 
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Fig. 2. Habitable zone lifetime as a function of stellar mass for different 
ambient stellar densities. The dotted line indicates the total lifetime of 
the host star, and the diamond indicates Earth. 



where ctt is the Stefan-Boltzman constant. Following Borucki 
et al. (2011), we calculate rJJz for an Earth-like planet using 
/ = 1, A = 0.3, and T cq = 307 K, which is the maximum 
temperature that allows for the presence of liquid water when 
accounting for the greenhouse effect. The luminosity is taken 
from stellar evolution models at solar metallicity (Marigo et al. 
2008). The results are shown in Figure 2 for ambient densities 
of E = {10 2 , 10 3 , 10 4 , 10 5 } pc~ 2 . The HZ lifetime decreases 
with density and stellar mass, due to the enhanced encounter rate 
and the large rg z , respectively. These estimates for the HZ life- 
time hold both for PPDs and planetary systems, since dynamical 
interactions would have comparable effects in both cases (e.g. 
Olczak et al. 2006; Parker & Quanz 2012). 

Figure 2 shows that, based on the Earth's existence alone, 
the solar system cannot have formed in a dense (E > 10 3 pc~ 2 ) 
environment, unless the ambient density decreased on a short 
(r is 100 Myr) timescale. Conversely, meteoritic evidence in- 
dicates that the young solar system must have endured nearby 
supernovae (e.g. Cameron & Truran 1977), which provides a 
lower limit to the product Et > 10 3 8 Myr pc~ 2 (see the re- 
view by Adams 2010). A plausible scenario is thus that the solar 
system formed in a massive (~ 10 4 M Q , E ~ 10 3 pc -2 ), but 
unbound association, which dispersed on a short (t ~ 10 Myr) 
timescale (see e.g. Dukes & Krumholz 2012, although they refer 
to such a system as a 'cluster'). Our results seem to disagree 
with Eisner & Carpenter (2006) who derive the disc fraction 
in the ONC and find no evidence of disc truncations. However, 
their conclusion may result from low-number statistics, and PPD 
mass estimates are more uncertain than radius measurements. 
On the other hand, our results agree with the studies of Bonnell 
et al. (2001) and Spurzem et al. (2009), and would explain why 
no planets have been found in the globular clusters 47 Tuc and 
NGC 6397 (Gilliland et al. 2000; Nascimbeni et al. 2012), where 
rE ~ 10 8 Myr pc~ 2 within the half-mass radius. This implies 
such a high number of encounters that it is improbable that any 
bound planets survived, most of them likely to have escaped the 
cluster due to two-body relaxation (e.g. Kruijssen 2009). 

We note that our theoretical estimates are conservative and 
provide upper limits to the disc sizes, because (1) we do not ac- 
count for potentially higher ambient densities in the past (e.g. 
Bastian et al. 2008), (2) we neglect the presence of massive stars 
at earlier ages of the SFRs, (3) we only consider the most dis- 
ruptive encounters and ignore the cumulative effect of weak per- 
turbations. Figure 2 thus provides strict upper limits to the HZ 
lifetimes. 
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Fig. 3. Evolution of the disc mass function (DMF) due to dynamical 
encounters at an age of 3 Myr for densities {2.25 x 10 4 , 2.5 x 10 3 , 5 x 
10 2 } pc~ 2 (solid, dashed, dash-dotted), sampling NGC 3603 at ra- 
dius intervals of {0" - 5", 16.5" - 43.5", 46.5" - 73.5"}. The dot- 
ted curve shows the initial DMF, and the vertical dotted line indicates 
the 4a sensitivity limit of ALMA Cycle 1 (30 min at 1.3 mm) at the 
distance of NGC 3603. 

4. Further observational avenues 

We present evidence for a change in the PPD radius distribution 
at ambient stellar densities of £ > 10 3 5 pc~ 2 at the > 97% 
confidence level, in line with the expected range due to close 
encounters with other stars on a ~ 1 Myr timescale. These den- 
sities are only reached in the densest parts of the ONC, which is 
not only consistent with the detection of reduced PPD masses in 
the centre of the region (Mann & Williams 2010), but also with 
studies concluding that encounters are not important in the ONC 
as a whole (e.g. Scally & Clarke 2001). Our results demonstrate 
that the stellar environment can be an important factor in setting 
the habitability of planetary systems. For instance, the existence 
of unbound, free-floating planets (see e.g. Bihain et al. 2009; 
Sumi et al. 2011; Strigari et al. 2012) is a natural outcome of 
our results. However, a ubiquity of Earth-like planets in the HZ 
of stars remains likely because a large fraction (~ 90%) of stars 
forms in unbound associations (see Kruijssen 2012 for a recent 
review, and observational references therein), of which the den- 
sity quickly decreases after their formation. 

To verify our results, more observations of PPDs in clus- 
tered environments are desirable. A fruitful approach would be 
to probe PPD truncations using disc mass measurements, which 
would provide much larger samples of PPDs. Figure 3 shows 
the effect of (only) dynamical encounters on the disc mass func- 
tion (DMF) using the disc mass loss description from Olczak 
et al. (2006) (see Appendix C). This will be easily observable 
in dense and young stellar clusters with ALMA. While the full 
ALMA array will be able to directly measure disc sizes, we pre- 
dict that, already from Early Science Cycle 1, the sensitivity will 
be sufficient to detect the variation of the DMF caused by the 
truncation (see Fig. 3), which would verify the result of Fig. 1. 
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Table A.l. Distribution of PPD sources over the host star type. 
'Unknown central star' indicates systems in which the central body 
has either not been detected or classified yet (mostly proplyd silhou- 
ettes in the ONC). Young stellar objects are classified as such if their 
enhanced envelope emission suggests a younger age than T-Tauri or 
Herbig Ae/Be. The 'KS' column shows the sample after applying den- 



Table A.2. Star-forming regions used in this study, listing the names of 
the regions, their numbers of objects N tj , and distances D. The fourth 
column shows the literature sources as (1) the Cores to Disks Spitzer 
survey (Evans et al. 2003), (2) the Taurus Spitzer survey (Rebull et al. 
2010), (3) the Robberto et al. (2010) survey of the ONC. 



sity and radius cuts for the KS test of Fig. 1. 




Name 




D(pc) 


Survey 








Lupus I 


20 


150 


1 


Type 


Sources 


KS 


Lupus III 


79 


150 


1 


Herbig Ae or Be 


15 


3 


Lupus IV 


12 


150 


1 


T-Tauri 


39 


21 


Ophiuchus 


297 


125 


1 


Unknown central star (ONC) 


36 


35 


Orion Nebula Cluster 


7759 


414 


3 


Young stellar object 


11 


8 


Perseus 


387 


250 


1 








Serpens 


262 


415 


1 








Taurus 


249 


137 
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Fig. A.l. Wavelengths A of the radius measurements for the final 
sample of protoplanetary discs versus ambient surface stellar density. 
Symbols have the same meaning as in the upper panel of Fig. 1. 



Appendix A: Properties of the sample 

After the selection procedure detailed in Sect. 2.1, we list the fi- 
nal PPD sample in Table A. 1 . To assess the heterogeneity of the 
sample, we show the wavelengths of the radius measurements as 
a function of ambient stellar density in Fig. A.l. The vast ma- 
jority of sources (75%) were measured in a narrow wavelength 
range below 3 /im (i.e near infrared wavelengths) with a spread 
of 0.25 dcx and centred at 1.32 /xm. The remaining sources were 
measured at millimetre wavelengths. For the ONC sample, 80% 
of all sources were measured at 0.66 /im with very little scatter 
overall. As shown in Fig. A.l, the 16 sources of the total sam- 
ple measured at millimetre wavelengths all have ambient surface 
stellar density < 200 pc~ 2 . Therefore the distribution at densi- 
ties above this value can be considered to be homogeneous. 

The regions from which the YSOs are taken to estimate the 
ambient density are summarized in Table A.2. Using these sam- 
ples, the ambient surface density of stars around each PPD is 
estimated as (Casertano & Hut 1985): 



S = 



N - 1 



(A.l) 



where N is the rank of the iVth nearest neighbour, and is the 
projected angular distance to that neighbour. We use N = 20, 
which is higher than the commonly-used value of N = 7 (cf. 
Bressert et al. 2010) and is chosen to improve the statistics of the 
density estimates. An additional effect of using a higher value of 
N is a slight decrease of the density estimates. This should be 
kept in mind when comparing our densities those in other work. 



Appendix B: A simple model for PPD truncations 

In this Appendix, we derive the upper limit to the radii of pro- 
toplanetary discs (PPDs) due to dynamical encounters. Where 
appropriate, we emphasize that the derivation is conservative, 
such that the obtained truncation radius is indeed an upper limit. 

Olczak et al. (2006) performed numerical simulations of 
disc perturbations and provided an expression for the relative 
disc mass loss AM/M due to encounters with other stars (their 
Eq. 4). If we assume that the mass loss occurs by stripping the 
outer disc layers and adopt the disc surface density profile of 
Sd oc r _1 used in their work, then Ar/r = AM/M. The 
expression for Ar/r from Olczak et al. (2006) is consistent to 
within a factor of three with the scenario in which a disc is al- 
ways truncated to the equipotential (Lagrangian) point between 
both stars. If the disc was already smaller than that radius, it is 
left relatively unperturbed. For the rough estimate made here, it 
thus suffices to write for the upper limit to the disc radius 



rd(r p ,m 1 ,m 2 ) 



\Jm 2 /mi + 1 



(B.l) 



where r p is the pericentre radius at which the perturber passes, 
mi is the mass of the perturbed system and m 2 is the mass of the 
perturber. The approximation of Eq. B.l follows Eq. 4 of Olczak 
et al. (2006) with reasonable accuracy for initial disc radii up to 
a few 10 3 AU (consistent with the parameter space in Fig. 1), en- 
counter distances r p > 0.002 pc (i.e. r p /rd > 0.2) and mass ra- 
tios W2/mi > 1. We have verified that these conditions are sat- 
isfied for the encounters that are expected to determine the disc 
truncation (see below). Following Binney & Tremaine (1987), 
the impact parameter b and the encounter radius due to gravita- 
tional focusing r v are related as 



2G(m 1 +m 2 ) 



v'r-n 



(B.2) 



where v is the relative velocity of the encounter. This equation is 
inverted to derive r p for each encounter. 

The truncation radius r^ of Eq. B.l depends on the variable 
set {b, v, mi, m 2 }, for which we specify probability distribution 
functions (PDFs). For the masses, we use a Salpeter (1955) type 
initial mass function in the range 0.1 M -m miH , where m max 
depends on age due to stellar evolution. For ages r < 4 Myr 
we assume m max = 100 M Q , while at later ages it is set by the 
Marigo et al. (2008) stellar evolution models at solar metallicity. 
The mass function is: 



dN 
oc - — oc m 
dm 



-2.35 



(B.3) 
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which is normalized such that J $ m dm = 1. Assuming a 
Maxwellian velocity distribution, the total number of encounters 
per unit velocity dv and unit impact parameter db follows from 
the encounter rate d 2 N/dbdv as (Binney & Tremaine 1987) 



oc Af(b,v) = t 



d 2 N 8n 2 vTb 



dbdv (47TCT 2 ) 3 / 2 



where v is the local number density of stars, r is the age of the 
region, and a is the velocity dispersion. The relative velocity 
ranges from v = O-oo and the impact parameter from 6 = 0- 
frfviAX (see below). As in Eq. B.3, we have normalized such 
that JJ $/vd&dv = 1, by writing $a/- = f dis Af(b,v)/n and 
defining n = JJ / d i s A/"(&, v)dbdv as the total number of encoun- 
ters at age r. The factor / d j s « 0.3 represents the fraction of en- 
counters that leads to disc mass loss according to Eq. B.l. This 
accounts for the fact that encounters with pericentres at inclina- 
tion angles 9 > 45° with respect to the disc plane cause only 
weak mass loss and retrograde encounters leave the disc almost 
unperturbed (Pfalzner et al. 2005). 

Given a sequence of encounters, the truncation of the 
PPD is set by the most disruptive encounter (Scally & 
Clarke 2001, although see Olczak et al. 2006), i.e. r Pi „ 



mm 7 ^min 



max 

). If we assume that {b 



are uncorrected, implying that the region is not mass- 
segregated, the PDF of the most disruptive encounter becomes 

^mini ^2, max .max ), (B.5) 

where P{b.v,m} represent the PDFs for the lowest b, lowest v and 
highest 77i2, respectively. Following the method of Maschberger 
& Clarke (2008, Eq. A5), these three PDFs are defined as 



/ /"bMAX 

Pb(bmin) = n<5> b (b min ) I J ^b(b')db' 



Pv(v m in) = n$ v (v min ) / § v (v)dv 



n-l 



(B.6) 



<E> m (m 2 )dra 2 



where &i,(b)& v (v) oc ^\r(b, v) are the distribution functions for 
b and v, with $b(6) oc vrb and § v (v) oc cxp (-v 2 /4a 2 )v 3 /er 3 , 
again normalized to unity in both cases. In Eqs. B.6, 6 m in, Vmin 
and m2 jmax indicate variable limits, and &max and m^MiN in- 
dicate fixed limits. The fixed limit buAX represents the max- 
imum impact parameter, which is given by the typical inter- 
stellar separation buAX — (48 /-kv) 1 ^ (the factor was chosen 
for consistency with Scally & Clarke 2001). It should be noted 
that while this is a physically motivated choice, it only weakly 
influences the result since the most likely most disruptive en- 
counter will typically be at 6 m j n -C &max- Assuming an age of 
t = 1 Myr, for surface densities of stars E < 10 5 pc -2 we find 
that pi, always peaks at impact parameters 6 m ; n > 0.002 pc (i.e. 
r p/ r d ~ 0.2), whereas for E > 10° pc~ 2 the most likely most 
disruptive encounter always has 777,2 > 0.5 M Q , which after av- 
eraging over the mass function to account for the distribution of 
mi gives 7772/mi > 2. This validates the use of the approxima- 
tion in Eq. B.l. 

By combining the Eqs. B.3 and B.5, the total PDF is 



$tot = Pb(&min)^(Wmin)jJm("72, m a X )$m(mi). 



(B.7) 



It should be noted that we did not include the mass of the per- 
turbed object mi in the PDF of the most likely most disruptive 



encounter (Eq. B.5), but instead average over the mass PDF it- 
self. The reason is that the stars in Fig. 1 span a range of masses, 
and a 'typical' relation between the truncation radius and ambi- 
ent density is preferable. 

Combining the previous equations gives a theoretical esti- 
mate for the typical truncation radius r tr as a function of the 
ambient density, velocity dispersion and age: 



I III™* 



1 dw mi nd777id7772. 1 



(B.8) 



where V indicates the complete phase space, i.e. 0.1 MQ-777 max 
in mass, 0-co in velocity and 0-&max in impact parameter. This 
expression provides the expected radius after the 'most likely 
most disruptive encounter', averaged over the stellar mass func- 
tion to account for the unknown mass of the perturbed system. 

Appendix C: Evolution of the disc mass function 

To calculate the evolution of the disc mass function (DMF), we 
assume that the initial disc mass 777d,i is related to the host stellar 

mass mi as 

m d ,i = /dmi, (C.l) 

where /d is a constant. We adopt /d = 0.03, which is in good 
agreement with observations (Andrews & Williams 2005) and 
sufficiently accurate for the order-of-magnitude estimate made 
in Sect. 4. Using a Salpeter (1955) stellar initial mass function 
(cf. Eq. B.3), the initial DMF is 

dJV(m d ,i) _ 1 dN( mi ) 

<3m d ,i /d dmi 

For each host stellar mass, we calculate the characteristics of the 
most likely most disruptive encounter as in Appendix B, using 
quantities that are appropriate for NGC 3603 (i.e. t — 3 Myr, 
g = 4.5 km s _1 , and R = 1.45 pc). Given a certain en- 
counter, the disc mass loss is calculated using the expression 
from Olczak et al. (2006, Eq. 4), which provides A = Amd/mj 
as a function of the host stellar mass mi, the mass of the per- 
turber 777 2 , the pericentre distance r p and the disc radius r&. 
To account for the dependence of A on the radius, it is calcu- 
lated for all radii from the observed sample at ambient densi- 
ties 10 _1 < E/pc~ 2 < 10 2 (see Fig. 1), including those with 
rd < 50 AU since the corresponding regions are all nearby and 
hence the detection limit is less stringent. At these densities the 
encounter rate is so low that the observed disc radii can be in- 
terpreted as 'initial' radii. The obtained values of A are then 
averaged to remove the dependence on r d , and integrated $ t ot 
(see Eq. B.7) in the same way as r tr in Eq. B.8. This provides 
the expected relative mass loss as a function of host stellar mass 
(A), 1 and hence the final disc mass is approximately 

777 d = (l-(A)) 777 d a. (C.3) 

The final DMF is then given by 

diV(777 d ) _ / d777 d \ _1 d7V(m dii ) 
d?77 d \dm d) i/ dm dj ; 

1 diV(777 d!i ) 



(1 - (A)) d777 d ,i 



(C4) 



1 Note that contrary to our Lagrangian approximation of Eq. B.l, the 
disc mass loss of the Olczak et al. (2006) equation does not increase 
monotonically with decreasing pericentre distance - for very close en- 
counters (typically r p /rd < 0.2) the disc mass loss is reduced. In such 
cases, the most disruptive encounter is not the closest encounter, and we 
account for this by adjusting r p to the value where (A) peaks. 
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